Localized spin ordering in Kondo lattice models 



O 

o 

> 
O 

o 



I. P. McCulloch^ , A. Juozapavicius^, A. Rosengren* and M. Gulacsi^ 
^Department of Theoretical Physics, Institute of Advanced Studies, 
The Australian National University, Canberra, ACT 0200, Australia 
^Department of Theoretical Physics, Royal Institute of Technology, 
SE-100 44 Stockholm, Sweden 
(February 1, 2008) 

Using a non-Abelian density matrix renormalization group method we determine the phase diagram 
of the Kondo lattice model in one dimension, by directly measuring the magnetization of the ground- 
state. This allowed us to discover a second ferromagnetic phase missed in previous approaches. The 
phase transitions are found to be continuous. The spin-spin correlation function is studied in detail, 
and we determine in which regions the large and small Fermi surfaces dominate. The importance 
of double-exchange ordering and its competition with Kondo singlet formation is emphasized in 
understanding the complexity of the model. 
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The Kondo lattice model (KLM) describes the inter- 
action between a conduction electron (CE) band and a 
half-filled narrow impurity, e.g. / electron, band and 
is thought to capture the essential physics of the rare 
earth compounds. Although intensively studied for two 
decades, the KLM is still far from being completely un- 
derstood. Recently, after the discovery of Kondo insula- 
tors and the non-Fermi liquid behavior, interest in this 
field has been greatly renewed, especially due to the non- 
Fermi liquid behavior discovered in most of the heav, 
fermion compounds, which resembles a Griffiths phase 

In order to understand the role of the impurity spin 
in determining the properties of KLM we must develop 
a better understanding of the magnetic correlations. 
The Griffiths phase in the one dimensional KLM oc- 
curs naturally;B it is, therefore, the prototypical model 
for heavy fermion compounds. Hence, this is an id| 
system to study since we have the bosonized solutio 
and we know the behavior of the CEs in both the param- 
agnetic (PM) and ferromagnetic (FM) phases. However, 
less attention has been given to understand the correla- 
tions between the impurity spins. This is the focus of our 
study. 

The Hamiltonian for the KLM is 



h = -* £ (4 



(i) 



where t > is the CE hopping parameter, Sj are 



The obtained result is m excellent agreement with a re- 
cent bosonized solutiono and contradicts common view 
that this phase boundary goes to infinite KondiQ|-jCOu- 
pling J as the CE density approaches half-filling.EfEl We 
also determine the regions of the phase diagram where 
large and small Fermi surfaces are dominant, which has 
been a central issue for much of the research in this area 
for some years. 

In addition, we have discovered a second FM region not 
seen before. For most dopings, this region of FM sepa- 
rates the regions of large and small Fermi surface. This 
most likely resolves the question as to the applicability of 
the Luttinger theorem to the KLM, shown by Yamanaka 
et a/.J9 since the Fermi points are not expected to remain 
constant across a phase transition. 

To accelerate the computation, we make use of sev- 
eral operators that commute with the Hamiltonian, 
5 + , 5~ , S z , I + ,I~ ,I Z , respectively the generators of the 
spin SU(2) and pseudospin SU(2) algebras.0 Combined, 
the generators form the algebra 50(4). All of the states 
in our DMRG calculation transform as irreducible rep- 
resentations of this algebra. Since 50(4) is non-Abelian 
these representations have, in general, degree > 1, which 
implies that a single basis state in the 50(4) repre- 
sentation is equivalent to multiple states of the purely 
Abelian representation of most previous DMRG calcu- 
lations. This is the origin of the dramatic performance 
improvements of the non-Abelian DMRG. The states are 
labeled by the eigenvalues of the Casimir operators of 



spin 1/2 operators for the localized spins, e.g. /, and 50(4), which are 5 2 = s(s + 1) and I 2 = 1). Hence 



= | J2a a' c j a (T <y-<y' c j a' with <x the Pauli spin matri- 
ces and Cj a , Cj g the electron annihilation and creation 
site operators. The Kondo coupling J is measured in 
units of the hopping t and partial conduction band fill- 
ing, n = N/L < 1, is assumed throughout. 

The method that we use is density matrix renormal- 
ization group (DMRG) which, however, is extended to 
explicitly preserve SU(2) spin and pseudospin symme- 
try. Hence we can measure the magnetization directly 
and determine rigorously the PM - FM phase boundary. 



we can label all irreducible representations by [s,i], which 
has degree (2s + 1) (2i + 1). In this construction, a chem- 
ical potential would appear as a term in the Hamiltonian 
proportional to I z , acting in an identical way a mag- 
netic field coupled to S z . Although the basis states in 
the calculation are eigenstates of 5 2 and I 2 , rather than 
S z and I z , all these operators mutually commute so it 
is possible to simply replace S z and I z by the chosen 
eigenvalues in this case. A single site of the Kondo lat- 
tice contains just three such states. The simplest is the 
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Kondo singlet state, transforming as the [0, 0] represen- 
tation of degree 1. The Kondo triplet state transforms 
as the [1,0] representation of degree 3, and encapsulates 

the three projections | frt), v^fl 1U> + I 4T», I U) 
in a single state. Here, ff denotes localized /, and J. 
the conduction electron spins, respectively. Finally, the 
holon state (actually, the tensor product of a holon and 
a / spin) transforms as the [1/2,1/2] representation of 
degree 4 and has the projections | fT 0), | ■!)■ 0), | fttl): 
| J|tl)- The single-site operators are 3x3 matrices over 
this basis. The matrix elements can be determined by 
the Wigner-Eckart theorem, which specifies the relation- 
ship between the 3 dimensional reduced basis and the full 
8 dimensional basis. For a comprehensive description of 
the new algorithm, see Ref. |[ At half filling (where the 
ground state is a pseudospin singlet) 400 block states are 
equivalent to around 2500 states of a calculation using 
N and 5 2 quantum numbers, although the relative ad- 
vantage of 50(4) decreases as the system is doped away 
from half filling. We used the new DMRG algorithm to 
obtain the ground state energy, magnetization and dif- 
ferent correlation functions, i.e., the momentum distribu- 
tion, density-density, conduction electron spin-spin and 
the / spin structure factor, S{k). The obtained results 
can be summarized with the phase diagram presented in 
Fig. 1, which will be analyzed in detail hereafter. The 
main properties of the phase diagram have been con- 
firmed on chains of 120 or more sites. Results for the 
magnetization were calculated on smaller chains, 40 - 60 
sites, where the energies can be calculated more accu- 
rately. We found no finite size effects that would affect 
the properties of Fig. 1. In all cases, we extrapolate to 
zero truncation error based on well-converged sweeps of 
between 200 and 500 50(4) states kept. 

As it can be seen from Fig. 1, the main feature domi- 
nating the KLM is / spin FM ordering. The FM ordering 
is due to the double-exchange (DE) interaction which ap- 
pears-as a consequence of an excess of localized spins over 
CEs:Elcach CEhas to screen more than one localized spin, 
and since hopping is energetically most favorable for CEs 
which preserve their spin, this tends to align the localized 
spins. This element was missing in the early approaches, 
which concentrated on the competition between Kondo 
singlet formation at large-/ and the RKKY interaction in 
the weak coupling limit £3 This picture is borrowed from 
the single impurity Kondo model and is inadequate for 
the lattice caseoB 

Starting the analysis of the phase diagram for large J, 
we see, that all CEs form singlets with the localized / 
spinsEil which become inert. The uncoupled / spins or- 
der FM in a mechanism similar to the J < case.El Here, 
there is no competition between Kondo singlet formation 
and DE. The fully polarized stateJwith 5 = [L - N)/2] 
appears for any value of n < lBO contrary to the sug- 
gestion of Refs. that close to half filling the PM phase 
extends to J — * oo. As J is lowered, KLM can be 
rigorously mapped into a random transverse field Ising 



model,0 hence the phase transition (the solid curve in 
Fig. 1) is identical to the quantum order - disorder tran- 
sition. It should be emphasized that this is also true for 
the second FM phase, as will be shown later on. 




FIG. 1. The obtained phase diagram of KLM. The two 
shaded areas are the FM phases. The open circles and tri- 
angles correspond to points at which the FM energy level 
crosses the S = level. The dashed curves are the derived 
phase transition lines (the solid curve was already obtained 
in Ref. |). 
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FIG. 2. Normalized magnetization curves (relative to the 
ground state energy, Eq) across the phase transition at quar- 
ter filling for a 40 site lattice. 

The phase transition obtained via DMRG fits excep- 
tionally well this picture, confirming the bosonization re- 
sult of Ref. |. The open circles correspond to points at 
which the energy of the FM state crosses the energy of 
the singlet state. Since the phase transition is second or- 
der, this is only an upper bound on the true transition 
line. However the partially polarized region is very small, 
of the order of J ft ~ 0.01, which is why this phase tran- 
sition has not previously been observed to be continuous. 
A typical example of the energy versus the magnetization 
(M) is presented in Fig. 2. This shows that in the tran- 
sition regime, d 2 E/dM 2 is positive. We have accounted 
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for all known random errors, these are errors arising from 
the tolerance of the matrix diagonalization, variations in 
the energy across the DMRG sweep, and error arising 
from the extrapolation to zero truncation error. These 
errors are of the order of the symbol size in this figure. 

Below the solid curve, Fig. 1, the Kondo singlets are 
not inert anymore and they greatly contribute to the 
properties of KLM. Excluding the Kondo triplet states, 
the CE wave function in the continuum limit satisfies a 
nonlinear Schrodinger equationtj which has finitely delo- 
calized solitonic solutions £3 This corresponds to a dress- 
ing of the CE by a cloud of antiparallel local spins, i.e., 
spin polarons are formed. The polaronic length scale 
competes with the length scale set by the free CE mean 
free path and introduces competing time scales: slow mo- 
tion of the polarons with low energy dynamics and fast 
motion of the free CEs with high energies. This sce- 
nario resembles a two-fluid picture with intrinsic inhomo- 
geneities which involves spin fluctuations and short-range 
spin correlations, which we call a polaronic liquid. 
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FIG. 3. Typical J dependence of the spin structure factor, 
S(k), and the momentum distribution, n(k) (n = 0.6). 

Finite temperature DMRG0 confirmed the presence of 
short-range / spin correlations in the van-Hove singulari- 
ties. Consequently the structure factor peaks at 2kp — n, 
where kp is the Fermi point determined by the filling of 
the CE band. This means that the localized / spins, even 
though they are completely immobile, contribute to the 
volume of the Fermi sea. This conventionally is called a 
Jarge Fermi surface, the effect of which is also seen in the 
momentum distribution function, see Fig. 3. As the po- 
larons are formed the peak of S(k) shifts from the small 
J/t value of 2kp: the slow motion of the spin polarons 
will dominate the low energy dynamics of the quasiparti- 
cles. This proves that the appearance of the large Fermi 
surface is a dynamical effect since it involves local inho- 
mogeneities, impurity spin fluctuation, and short-range 
correlations of the / spins. This phase is relatedo to a 
Griffiths phase, suggesting that the small - large Fermi 
surface crossover is a Griffiths singularity. 

The large Fermi surface is conventionally explained 



by reference-to the periodic Anderson model (PAM) 
ancestryBErE-3 Our results imply that even for PAM, this 
simple picture is inadequate. In particular, we see no rea- 
son why a small - large Fermi surface crossover, marked 
by a FM phase, should not also appear in PAM. How- 
ever, the behavior of the Fermi surface crossover close to 
quarter filling is numerically difficult to determine (dot- 
ted line in Fig. 1), hence we are not yet able to rule 
out the possibility that the large and small Fermi surface 
regions are adiabatically connected. Even prior to the 
current calculation, the nature of the Fermi surface in 
the weak-coupling regime was not clear, with the sugges- 
tion from Ref.U that the Fermi surface vanishes at a point 
in proximity to where we find the ferromagnetic phase. 
For n < 0.5 the width of the polarops is over several lat- 
tice spacings (diverging for n — > 0|I3) hence the energy 
needed to excite these polarons is too large for this effect 
to happen. The polarons will not contribute to the low 
energy dynamics and the system behaves as an RKKY 
liquid, as we explain below. 

An interesting phenomenon appears as we further 
lower J. The residual weight attached to the Kondo sin- 
glets vanishes, hence all CEs which participated in the 
formation of these singlets, become delocalized. The dis- 
tance between these CEs is much larger than the lattice 
spacing, and below J < 2y / nsin(7rn) their continuum 
limit takes the regular quantum sine-Gordon form.0 In 
the bosonization language of Ref. |[ this means that the 
spin Bose fields, $o- cannot be approximated by their 
noninteracting expectation values, rather by their expec- 
tation value corresponding to a sine-Gordon (sG) model, 
(f> CT w (^ a ) S Q. However, the charge degrees of freedom 
not being affected by the sine-Gordon spin gap, their 
corresponding Bose fields, $ p may be still approximated 
by their noninteracting values. Extending the bosonized 
results of Ref. || to a finite ($ ct ) s g, we obtain the crit- 
ical Hamiltonian governing the PM - FM phase tran- 
sition at intermediate J values in the form: -ff cr jt = 
-J 2 ^1/(2^V)E j S|S| +1 +2J6E j {1-((^> s 2 g/2)[1 + 
J/(2itvf)] 2 + cos(2kpj)}S : j , where A and B are func- 
tions which depend only.-nii the cutoffs introduced by 
the bosonization scheme, H Following closely previous 
bosonization approaches, Bp we can prove that the critical 
behavior of the FM transition for the intermediate this 
J case is of a random transverse-field Ising model type, 
where the transverse field hj = 2JB{1 - (($ CT )g G /2)[l + 
J/(2itvf)] 2 + cos(2kpj)} is driven by a displaced co- 
sine distribution of the form: p(h) = [1/ (2irJB)]{l — 
[h/(2JB) + «*«r)k/2)[l + J/{^v F )Y ~ I] 2 }- 1 ' 2 . Ac- 
cordingly, the FM transitions emerging at intermediate 
values of J are of a quantum order - disorder type. 
These transitions are driven by spin polarons, contrary 
to the FM phase emerging at high J values, which is 
given by the uncoupled / spins (in a mechanism simi- 
lar to the J < case). The new critical line is J c = 
a{A, B) sin(7rn/2)/[l - P(A, B)] - j{A, B, ($ CT ) S 2 G ). The 
bosonization (conformal field-theory) arguments docs not 
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determine the magnitude of a, and 7, accordingly these 
constants are used as fitting parameters to the numeri- 
cally obtained points. The best fits are the dashed curves 
in Fig. 1. 

This is the second FM phase in Fig. 1, which has 
proven .difficult to detect with conventional (Abelian) 
DMRG-llaO Previous DMRG calculations did show a 
weak FM signal at n = 0.8 and J = 1.6 and 1.8,E3 but 
the results were discarded in later papers by the same 
authors .Era Likewise an exact diagonalization of a vecy. 
small system gave FM for n = 0.75 and J = 1.5.E3 
Using the non- Abelian DMRG algorithm we could also 
check the energy of each total spin state, shown in Fig. 
4, which clearly shows a second ferromagnetic region al- 
though we have not yet confirmed numerically the or- 
der of the phase boundaries. For the FM Kondo lattice 
model, J < 0, a phase separated regime was observed 
in numerical approaches .113 However for J > we found 
no change in sign of the inverse charge compressibility. 
Thus, this phase is a true FM rather than a phase sepa- 
rated region. 
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FIG. 4. The gap, AE, from the fully polarized ferromag- 
netic state to every other spin state vs. J, for n — 0.8, and 
60 site lattice, measured along intervals of J ±0.05. For most 
data points the error bars are of order ctab ~ 10" J or less, 
except for the S — curve for very low and very high J, where 



the errors are of order ctae r 
second ferromagnetic region. 
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Below the second FM region the KLM reduces to a 
system of free localized spins in fields determined by CE 
scattering: dominant 2kp modulations are manifest, see 
Fig. 3, superimposed on an incoherent background. This 
reflects the momentum transferred from the CE band to 
the spin chain in backscattering interactions, together 
with incoherent forward scattering. This case is referred 
to as an RKKY liquid as the scattering processes give an 
RKKY-like correlation for the / spins, even though the 
RKKY interaction strictly diverges in one dimension. 

In conclusion, using a non- Abelian DMRG method a 
most comprehensive analysis of the short- and long-range 
ordering of the localized moments in KLM is presented. 
We show that DE ordering and its competition with 



Kondo singlet formation is the dominant feature of the 
phase diagram. The non- Abelian DMRG method allowed 
us to discover that FM does not only appear at large J 
but also at intermediate values. This second FM phase 
was missed in previous approaches. We also show that 
at large J FM is due to ordering of uncoupled / spins, 
while for intermediate J, i.e., the second FM region, FM 
is due to ordering of the spin polarons. The inhomoge- 
neous polaronic state between these two FM phases is 
analogous to a two-fluid system and it exhibits a large 
Fermi surface. 
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